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Abstract 

NGC3718 is a LINER L 1.9 galaxy, lying at a distance of about ~ 17.4 Mpc away from earth and its similarities with NGC5128 often 
award it the name “northern Centaurus A”. The presence of a compact radio source with a candidate jet structure, a prominent dust 
lane, and a strongly warped molecular and atomic gas disk are indicative that NGC3718 has undergone some sort of a large scale 
gravitational interaction sometime in the recent past. This channeled gas towards the center, feeding the black hole and igniting the 
central engine. One proposed scenario involves an encounter with the close neighboring galaxy NGC3729, while other authors favor 
a merging event with mass ratio > (3 - 4) : 1, as the origin of NGC3718. 

We use high angular resolution (~ 100 mas) e-Merlin radio and SUBARU NIR (~ 170 mas) data, to take a detailed view of the 
processes taking place in its central region. In order to preserve some objectivity in our interpretation, we combine our results with 
literature values and findings from previous studies. Our NIR maps suggest, on one hand, that towards the stellar bulge there are no 
large scale absorption phenomena caused by the apparent dust lane and, on the other, that there is a significant (local) contribution 
from hot (~ 1000 K) dust to the nuclear NIR emission. The position where this takes place appears to be closer to the offset compact 
radio emission from our e-Merlin 6 cm map, lying offset by ~ 4.25 pc from the center of the underlying stellar bulge. The shape 
of the radio map suggests the presence of one (or possibly two, forming an X-shape) bipolar structure(s) ~ 1 (~ 0.6) arcsec across, 
which combined with the balance between the gas and the stellar velocity dispersions and the presence of hard X-ray emission, point 
towards effects expected by AGN feedback. We also argue that NGC3718 has a “core” in its surface brightness profile, despite the 
fact that it is a gas-rich galaxy and we discuss its mixed photometric and spectroscopic characteristics. The latter combined with the 
observed spatial and radio offsets, the relative redshift between the broad and the narrow Ha line, the limited star formation activity 
and AGN feedback, strongly imply the existence of an SMBH recoil. Finally, we discuss a possible interpretation, that could naturally 
incorporate all these findings into one physically consistent picture. 

Key words. Galaxies: kinematics and dynamics - Galaxies: active - Galaxies: evolution - Galaxies: formation - Galaxies: photometry 


1. Introduction 

Knowledge of the circumnuclear activity in active galaxies is es¬ 
sential for understanding the fueling of the central engine. Our 
understanding of the underlying physical processes that con¬ 
tribute to the nuclear activity is, however, still far from being 
complete and conclusive. While on large scales (> ~ 3 kpc) 
the picture is clearer, as large scale dynamical perturbations i.e. 
galaxy collisions, mergers (e.g. Toomre & Toomre fl972| ) etc., 
are proposed as the mechanisms responsible for removing an¬ 
gular momentum from the gas, driving it towards the central re¬ 
gion, the corresponding processes at smaller scales (sub-kpc) are 
not very well understood. Mechanisms ranging from nested bars 
(e.g. Shlosma n et al.| 1989) and spirals (e.g. Martini & Pogge 
1999) to warped nu clear disks (Schinne rer et al.| ; 2000b a) and 
m - 1 instabilities ( Kormendy & B ender 1999; Garcia-Burillo 


et al.|2000| , have been proposed over the years in order to ex¬ 
plain the smaller scale phenomena, the discussion, however, is 
still open. 

An ideal laboratory for trying to shed light on some of 
these important questions is NGC3718 (figure[l]). NGC3718 and 
its (supposed) companion NGC3729 belong to the loose Ursa 
Major group, with NGC3718 being one of the largest galaxies 
in the group. It is unclear, however, whether these two galaxies 
intera c t gravitationally or not and to what extent ( [Tully et aL 
|1996| ). [Karachentsev et al7j |2013| find that both galaxies be¬ 
long to the sub-group NGC3992 (named after the gravitation¬ 
ally dominant member), which is one of the most massive sub¬ 
groups of Ursa Major. They also note that, velocities within the 
NGC3992 sub-group do not show any visible correlation with 
distances, indicating a non-virialized system. Morphologically, 
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NGC3718 is classified as an SB(s)a pec by de Vaucouleurs et al 
( |1991| ), mainly due to the prominent -’’spiral arm” ending- dust 
lane which was consid ered to be indicative of the presence of 
a bar. |Holmberg| ( |1958[) classified NGC3718 as an SOp, whereas 


other authors (e.g. Whitmo re et al.|19 90; Reshetnikov & Combes 
|1994| ), classify it as a polar ring galaxy, confusing our under¬ 
standing of its true morphology. 

The distinguishing features of NGC3718 are the prominent 
dust lane, which runs across the entire stellar bulge, and its 
strongly warped molecular and atomic gas disk. Several au¬ 
thors have extensively studied the gas dynamics of NGC3718. 
Schwarz ( |1985| ) studied the HI dynamics and found that the 
atomic gas distribution forms a 3-dimensional warped structure, 
which could be described by tilted but concentric rings, orbit¬ 
ing from n early edge-on at sma ller radii, to nearly face-on at 
larger radii. |Pott et al.| ( [2004] ) and |Krips et al.| ( [2QQ5 ) studied the 
molecular gas distribution using CO( 1 —> 0), CO( 2 —> 1) and 
HCN( 1 —> 0) as tracers. They successfully fit tilted rings on 
NGC3718 and found that the molecular gas motion, generally, 
follows that of the HI gas but the warp continues down to scales 
of ~ 250 pc. Sparke et al. ( ]2009| ) re-mapped the HI distribution 
using higher resolution VLA data and confirmed (though with 
slightly different parameters) the tilted ring models of the afore¬ 
mentioned studies. They show that the outer gas orbits extend to 
~ 35 - 42 kpc and they estimate its age at ~ 2 - 3 Gyr, whereas 
the inner gas orbits are nearly polar and still under formation. 
Finally, they do not see any HI gas in the plane of rotation of 
the stellar disk and they argue in favor of the classification of 
NGC3718 as a polar ring galaxy, invoking differential preces¬ 
sion in order to account for the warp. 

Some of these features are also present in NGC5128, host 
of the famous Cen taurus A radio source. A simila r dust lane 
([Dufour et al.|1979] ) and a warp in the gas disk (e.g. Wild et al. 
|1997[|Spa rke 1996 ), along with the presence of radio emission 
(though more dominant) from the nucleus, are common proper¬ 
ties which often lead to the characterization of NGC3718 as “the 
Northern Centaurus A”. 

NGC3718 is also one of the NUGA sources, a survey aimed 
at the study of nearby low-luminosit y active galactic nuclei 
(LLAGN) ( Garcia-Burillo et aL|2003| . Spect roscop ically, it is 
classified as a LINER type LI.9 galaxy ( |Ho et al.]|1997b| ). A 
weak broad Ha emission component with FWHM Ha ^road ~ 
2350 km s -1 is detected, originating from the nucleus. 
Additionally, the presence of strong [0/]/L16300A with 
FWHM\ 


m 


570 km s 1 is indicative of a hidden AGN 


( [Filippenko & Sargent|[1985| ). |Krips et af| ( j2007[ ) detect a can¬ 
didate jet structure in a 18 cm Merlin radio map, lying NW of 
the nucleus and stretching to ~ 0.5 arcsec. They also measure 
the bolometric luminosity of ~ 10 41 ergs s -1 which implies a 
sub-Eddington system. 

|Chitre & Jog| ( |2002| ) treat ARP214 (an alternativ e name 
for NGC3718) as an advanced merger remnant, while Jog &| 
|Chitre| ( [2002| discuss its mixed characteristics. Photometrically, 
it shows an exponential light profile (like a spiral galaxy) but 
kinematically, it is mainly supported by pressure from random 
motion of stars (like an elliptical galax y), as indicated by the 


V C 


0.5-1 within the inner few kpc. 


Bournaud et al. 


(2005) 

lane- 


simulated galaxy mergers with mixed photometric and 
matic characteristics through N-body simulations and they find 
that such objects could result from mergers with mass ratios 
>(3-4.5): 1. 

This paper is organized as follows: In section [2] we present 
the data and data reduction, in section [3] we describe the pre¬ 


processing and alignment procedures, in section[4]we present the 
NIR maps, we attempt a light decomposition and we present the 
6 cm e-Merlin radio map, in section[5]we discuss the scaling re¬ 
lations, the classification and the observed mixed characteristics 
of NGC3718, in section [6] we present additional observational 
evidence regarding the presence or not of an SMBH recoil and 
in section [ 7 ] we try to put all the pieces together in order to get a 
physically consistent picture. Finally, in appendices |A| [Bl and |C| 
we provide additional information for the various arguments we 
present in this paper. 

Throughout this paper we adopt the following: In all im¬ 
ages North is up and East is left. The cosmology values used 
are Cl M = 0.21, Q A = 0.73, H 0 = 67.8 km s -1 , with a redshift 
of z = 0.003927 corrected to the reference frame defined by the 
3K CME0 These suggest a distance modulus fi - 31.2 and a 
distance D = 17.4 Mpc for NGC3718. The cosmological scale 
for the adopted cosmology at this distance is 84 ——. We also 
use the “Mass to Light” ratio relationship from Bell et al. (2003) 
for all the conversions from light to mass. The general form of 
their j- ratio in the NIR and optical (B-V) colors is: 


log 


MJM Q \ 

l/l^ a I 


- a A + b A (B - V). 


( 1 ) 


with a A and b A being aj = -0.261, a H = -0.209, a K = 
-0.206 ,bj = 0.433, b H = 0.210, and b K = 0.135. For NGC3718 
this value is ^ ~ 0.78, using My GC37ls = -20.73 mag and 
MNGC 3718 _ _20.01 mag, magnitudes corrected for extinction 
and K-corrections, taken from NED (footnote [I]). 


2. Observation and data reduction 

2.1. SUBARU data 


The data set used for the analysis of NGC3718 consists of, AO 
assisted near infrared (NIR) data in J (5), H (11) and K s (12) 
bands (number of images) respectively, taken on the 17th of May 
2012, with the SUBARU telescope at Mauna Kea, Hawaii, using 
the HiCIAO ( [Suzuki et al. 2010| ) instrument, operating in the 
Direct Imaging Mode. The individual frame exposure time for 
all bands is t J e ^ s = 60 sec, with the total duration of the ob¬ 
servations being ~ 1 h 40 m. During this time, the total variation 
of airmass is ~ 0.082. The AO 188 AO system is used (Hayano 
et al.||2010l ). It is equipped with a 188-element wavefront cur¬ 


vature sensor with photon counting APD modules and a 188 el¬ 
ement bimorph mirror, installed at the IR Nasmyth platform of 
the Subaru telescope, which for this dataset operates in self ref¬ 
erence mode on the core of NGC3718. As a result, the angular 
resolution of the data is ~ 170 mas. The 2048 x 2048 pixels 2 
Hawaii-IIRG HgCdTe detector provides a pixel scale of 0.010 
ar p c l e e i , with a FoV (Field of view) of 20 x 20 arcsec 2 . 

No reduction package was available for HiCIAO, so a 
pipeline was developed from scratch, in order to correct for the 
high frequency 32-strip artifact noise, introduced by the 32 read¬ 
out channels of the detector. All images were dithered and have 
undergone the usual bad pixel correction, flat-fielding (dome- 
flat), alignment (based on ellipse fitting) and median coadd treat¬ 
ment. The final images allow us to generate sky frames for 
each band from the data themselves, by clipping them near the 
modal background value, and adding a high frequency layer on 


1 http://ned.ipac.caltech.edu/cgi-bin/objsearch? 
search_type=0bj_id&objid=2688Q&objname=l&img_stamp= 
YES&hconst=67.8&omegam=0.27&omegav=0.73&corr_z=l 
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Figure 1. 9.5 x 4.5 arcmin SDSS gri composite image of 
NGC3718. 


to flfig ure [ 2 ]). The data with the sky subtracted are shown in 
figure [5] 



Figure 2. K s band before sky subtraction (left) and the sky frame 
extracted from the K s image (right). Similar sky frames were 
created and subtracted from the J and H bands, as well. 


2.2. e-Merlin data 

The e-MERLIN synthesis telescope is a seven element inter¬ 
ferometer, with baselines of up to 217 km and connected by 

2 Perhaps a more correct term would be a “background” subtraction, 
but in the NIR the majority of the background contamination refers to 
the sky component. We create our “sky” frame(s) by keeping values 
lower than the mode of the image(s), which should be a good approxi¬ 
mation of the sky contamination. On top of that, we add a layer obtained 
by applying a high pass filter on the image(s), in order to account for 
the pixel-to-pixel variations as well. 


J Band h Band 



Figure 3. /(top left), //(top right) and K s (bottom middle) ready- 
for-science images. 


a new optical fiber network to Jodrell Bank Observatory near 
Manchester, UK. An inhomogeneous array, e-MERLIN is com¬ 
prised of the 76 m Lovell telescope, a 32 m dish at Cambridge, 
and the following 25 m antennas: Mark II, Knockin, Defford, 
Pickmere and Darnhall. 

The data in this work were taken during the commissioning 
phase of e-MERLIN with only the five 25 m dishes, resulting in a 
primary beam of 10', a maximum baseline of 133.7 km (between 
Pickmere and Defford) and a minimum baseline of 11.2 km (be¬ 
tween Pickmere and Mark II). The final, fully expanded, array 
will have a bandwidth of 2 GHz, providing more than 10 times 
the continuum sensitivity as the original MERLIN. NGC3718 
was observed at 5 GHz for 10 hrs by e-MERLIN on the 3rd of 
August 2011, with 4 128 MHz sub-bands of 512 spectral chan¬ 
nels each, yielding a total bandwidth of 512 MHz. The observa¬ 
tions did not include baselines to Cambridge. The final angular 
resolution we obtain is ~ 100 mas. 

Data were reduced and analyzed using the National 
Radio Astronomy Observatory’s (NRAO) Astronomical Image 
Processing System (AIPS). Data were initially edited with 
SPFLG and IBLED, averaged to 64 channels per intermediate 
frequency (IF) channel and concatenated with DBCON before 
further editing was conducted. FRING was used to derive de¬ 
lay and rate corrections for the calibrator sources, while CALIB 
was used to derive time-dependent phase and then amplitude and 
phase solutions. Flux calibration was performed using short ob¬ 
servations of 3C286 at the beginning and end of each run and the 
flux density scale was calculated with SETJY ( [Perley & Butler 


2013). The flux density for each IF was then reduced by 4%, 


in order to account for the resolution of the e-MERLIN shortest 
spacing (see the MERLIN User Guide). Bandpass and phase cal¬ 
ibration was performed using the bright point sources. Absolute 
calibration is expected to be accurate to about ~ 10% for com¬ 
missioning data. 

The quoted positional errors, associated with phase trans¬ 
fer errors from phase reference to target source (for a typical 
3° separation), are ~1 mas. This assumes 6-stations including 
Cambridge and online L-Band link corrections. The phase ref- 
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erence source used in our e-Merlin observations is J1146+5356 
and we have referenced our observations to the coordinates given 
in the VLBA calibrator list, namely tr 2 ooo=11:46:44.204328 and 
£> 20 oo=+ 53 : 56 : 43 . 08356 . This places J1146+5356 about 2.28° 
away from NGC3718. Taking into account the lower resolution 
of these data without Cambridge, and making some allowance 
for known small phase corrections for the L-Band link trans¬ 
mission delays around our observational epoch, we conserva¬ 
tively estimate the phase transfer errors to be around ~4 mas. 
The phase reference source position from the VLBA Calibrator 
list quoted formal positional uncertainty errors of 0.35 mas 
and 1.17 mas in R.A. and Dec., respectively. The formal mea¬ 
surement error from the peak of the e-Merlin image will be 
100/2Q mas with a restoring beam of 100 mas and a peak signal 
to noise in the image of Q which is better than 60 in our case. 
Hence the total estimate of the positional uncertainties result in 
a value of ~4.1 mas. 


3. Data processing 

3.1. Centering 

The first step of our SUBARU data analysis, is to align the J,H 
and Ks images. In the absence of secondary sources in our FoV, 
we choose an alignment based on ellipses-on-isophote contours 
fitting. For this purpose we use IRAF’^Jtask ellipse , which fits 
ellipses to iso-intensity contours of a galaxy’s light distribution. 
We plot the X center and Y center coordinates indicated by each fitted 
ellipse versus its S(emi) M(ajor) A(xis), as shown in figure [4] 


tively normall}j^] The behavior of the J band, can be attributed 
to the lower, with respect to the H and Ks bands, total flux, Fj 
(. F Ks ~ 3 Fj and F H ~ 2.5 Fj respectively). This is partly due 
to the lower integration time (tj = 300 sec, t H = 660 sec and 
t Ks = 720 sec), but it mainly indicates that the stellar flux is 
of the order of the one expected, from an evolved stellar popu¬ 
lation of K/M type star^] The lower flux results in more noisy 
isophote contours and, consequently, in higher uncertainties. The 
dust lane of NGC3718 does not seem to affect the light distri¬ 
bution in the NIR to a large extent, mainly because it lies far 
away enough (see section [4] and Appendix [A|. Moreover, the 
dominant component, both in terms of structural size and illu¬ 
mination, should be the stellar bulge of the galaxy. The absence 
of lame scale contour deformations, and the general contours 
shapqj indicate a largely relaxed stellar component in all three 
NIR bands. 

We therefore consider, that the average of the most promi¬ 
nent stable values of the (. X center , Y center ) coordinate pairs (points 
within the shaded ellipses of figure]?]), should act as a very good 
first order approximation of the photocenter (i.e. center of the 
stellar bulge) for each band. The uncertainty of each (X center , 
Ycenter ) pair, is considered to be the quadratic addition of the stan¬ 
dard deviation of the most prominent stable values used for each 
band with their average fitting uncertainties ( dX center , dY center ) 
as computed by the ellipse task itself. The suggested photocen¬ 
ters and their corresponded uncertainties, are displayed in table 

m 

Table 1. First order approximation photocenters. 



SMA (pixel) SMA (pixel) 


Figure 4. X (left column) and Y (right column) central coordi¬ 
nate of each fitted ellipse versus its SMA, for the K s (red/top), 
//(blue/middle) and /(green/bottom) bands. The shaded ellipses 
indicate the most prominent stable values used to derive the first 
order approximation photocenters. 


Apart from the X center coordinate of the J band and the Y center 
coordinate of the Ks band, the rest of the curves behave rela- 


3 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 


Filter 

y 

^center 

(Pixel) 

y 

± center 

(Pixel) 

r v total 

center 

(mas) 

T -rrtOtCll 

a * center 

(mas) 

J 

903.16 

1058.36 

3.4 

2.2 

H 

897.10 

1062.59 

2.5 

3.0 

K s 

882.39 

1059.33 

1.5 

1.3 


Notes. First order estimation of X center & Y center and the uncertainties 
involved in mas, as suggested by ellipse fitting on the 7, H and K s 
bands. 1 pixel equals 10 mas. 


3.2. Pivot and Subtract 

Having a first order approximation of the center of the stellar 
bulge, the next question is: “What could be the cause of the devi¬ 
ation in the central coordinates of the K s band?”. A clean bulge’s 
light distribution should have a nice symmetrical bell shape 
structure (e.g. a Sersic profile), as are all spherical/elliptical 
shapes that are projected on a 2-D surface. This should lead 
to co-centered isophotes and, therefore, to co-centered fitted el¬ 
lipses. 

So we consider the following scheme: A theoretical perfectly 
symmetrical 3-D structure should look the same under any rota- 

4 In the absence of asymmetries in the light distribution, the plots of 
figure [4] are expected to be close to horizontal lines i.e. all the fitted 
ellipses should point towards a common center. 

5 The V - K color of AGC3718 is ~ 3 (NED), suggesting a dominant 
stellar population with T eff ~ 4000 K, consistent with K/M type stars 
(Ridgway et al.|1980j). 

b The NIR ellipticities range between ~ 0.09 - 0.12 in J band, ~ 
0.07 - 0.09 in H band and ~ 0.03 - 0.07 in Ks band, in the radial 
interval between ~ 0.5 - 1.5 arcsec. 


4 

















Markakis et al.: NGC3718 - inside the heart of a LLAGN 


tion around the perpendicular to the X-Y plane axis, if we have 
a spherical shape, or at least under n x n rotation, with n being 
an integer, if we have an elliptical shape. So if one was to sub¬ 
tract a 2-D projection of such structure (a bell shape curve) from 
a rotated (around its center) copy of itself, one should receive 
0 as a result. Following this scheme for a perfectly symmetri¬ 
cal, isolated and undistorted bulge, if one rotates and subtracts 
its 2-D projection from itself, one should not receive any ob¬ 
vious residual pattern, apart from random noise. In the case of 
NGC3718 we have an elliptically projected bulge, so this scheme 
is valid only under nxn rotation around the first order approx¬ 
imation photocenters, derived from the curves of figure [4j with 
n- 1,2, 3, ...n. 

The pivot (n rotation) and subtract operation reveals strong 
residuals, especially in the K s and H bands (figure [5]), whereas 
in the J band the residuals are considerably more noisy. 


I Band H Band 



Figure 5. /(top left), //(top right) and (bottom middle), 
smoothed with a 10 pixel Gaussian kernel, residual images from 
the pivot and subtract operation. The residual mean peak fluxes, 
expressed as a percentage of NGC31 18’s mean peak fluxes, are 
~ 4.5%(/), ~ 3.0% ( H ) and ~ 3.4% (K s ) respectively. 


3.3. Evaluating the centering 

Despite the fact that the uncertainties implied are of the order 
of a few mas (see table p}, we test the robustness of our initial 
center estimations, by shifting the pivot point in all bands by 
± 1. Y^ter ± D (.X;± 0-5, Y^ ter ±0.5 ) and repeat- 

ing the pivot and subtract operation. In this way, we test a series 
of sixteen alternative pivot points, arranged in a rectangular (2 
pixels in side) pattern around the initial center estimations. If the 
mean central (300 x 300 pixels 2 ) residual produced by a can¬ 
didate (X 1 ter , ycenter ) * s smaller than the mean residual of the 
(XCenter , ycenter)? then it is adopted as a better center estimation 
and the process is iterated, until the minimum mean residual is 
reached. 

Our initial center estimations proved to be quite accurate, at 
least within ±0.5 pixels (±5 mas) per coordinate, since these are 


the pivot points that produce the minimum mean residuals. This 
is slightly larger than the statistical uncertainties of table [T] and 
we adopt them, as a more conservative approach. An example of 
the above process for the H band, is shown in figure[6] In general, 
the residuals produced by the alternative pivot points are a lot 
less well defined, especially in the central region. The brighter 
central features that appear are indicative of the, deliberately, 
incorrect pivot points used. The extended light excess towards 
the NE from the center, is largely unaffected, so the larger scale 
asymmetries can still be traced and removed. 



Figure 6. Examples of the H band smoothed residual images 
after a rotation around incorrect pivot points, by (X l " e lt nter - 
0.5, Y^ ter + 0.5) (upper left), (X^ ter + 0.5, Y^ ter + 0.5) (upper 
right), (X^Lr Y ^ter) (center), ( X^ ter - 0.5, Y^ ter - 0.5) (lower 
left) and (X l c n e * ter + 0.5, Y l ^ nter - 0.5) (lower right). 


We remove the residuals of figure[5]from the images of figure 
[3] in order to clear the light distribution from all the asymmetries, 
as shown in figure [7] The symmetrical light images in all bands 
are clearly more round and generally, better defined. 

We note that the residual images represent asymmetries be¬ 
tween diametrically opposite parts of the light distribution. In the 
presence of inclination and/or tilt, the further we move from the 
center, the line-of-sight light path (and therefore the brightness) 
difference between diametrically opposite parts of the galaxy be¬ 
come increasingly larger, so that, any large scale asymmetries 
can not be simply treated as light over-densities with a physical 
meaning. As we approach the center however, these effects be¬ 
come increasingly smaller, so that, in small scales the observed 
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asymmetries should mainly represent true light (and therefore 
mass) over-densitie^] The latter suggests that we can safely use 
the residual maps and the symmetrical light images for extract¬ 
ing physically robust quantities, as long as we restrict ourselves 
in small scaled 

J Band H Band 


Ks Band 



Figure 7. /(top left), H( top right) and K s (bottom middle) sym¬ 
metrical light images after the removal of the residuals. 


3.4. Angular resolution matching & flux calibration 

Having an accurate center estimation, we align our images with 
respect to the / band. We match the angular resolution of our 
images, by convolving the better seeings (H and K s ) to the 
worst one (/), using the FWHMs as measured from our PSF- 
reference stars in each band. The Gaussian kernel used, equals 
to the quadratic difference of the smaller from the larger FWHM, 
as follows. 

S F whm(H) = V[FWHM(J) 2 - FWHM(H) 2 ] * 9.8 (2) 

6fwhm(K s ) = VtFWHM(J) 2 - FWHM(Ks) 2 ] ~ 10.3 (3) 

Flux calibration was performed, using our standard stars, to de¬ 
rive Zero Points for each band according to: 

mi = ZPi - 2.5 log([countsi/Exp.Timei]) (4) 

3.5. Coordinate uncertainties 

In order to translate the relative into an absolute coordinate sys¬ 
tem, the following method is used: We align our SUBARU / 

7 As we approach the center, the brightness difference becomes in¬ 
creasingly smaller. When the residuals are minimal, with respect to 
neighboring pivot points, they should represent, to a first order, actual 
light density variations. 

8 Following Ockham’s razor, however, throughout this paper we try 
to restrict the use of either of these to minimum, for the sake of credi¬ 
bility. 


band with a re-sampled and re-scalecQ version of the SDSS z 
band image, in order to (astrometrically) calibrate it. We choose 
the SDSS z band, because its central wavelength (0.9134/im) is 
the nearest to our / band (1.220/mi), so the centroids of similar, 
old stellar populations are expected to coincide in the NIR, (see 
Appendix [B] and footnote [5]), and also because the SDSS is the 
most_ accurate (in astrometric terms) survey currently available 
( |Pier et aL||2003| . The total intrinsic SDSS astrometric uncer¬ 
tainties were calculated following the instructions described in 
the SDSS online documentation. The resulted values for our en¬ 
tire uncertainty budget, are shown in table [2] 


Table 2. Astrometric uncertainties and coordinates. 


uncertainties: SDSS z-J 

Cent./Shift e-Merlin 

Total 


(mas) (mas) 

(mas) 

(mas) 

(mas) 

J band 

6(R.A.) 

42.5 10.8 

5.0 

4.1 

44.3 

6(Dec.) 

41 11.0 

5.0 

4.1 

42.9 

H band 

6(R.A.) 

42.5 10.8 

7.1 

4.1 

44.6 

6(Dec.) 

41 11.0 

7.1 

4.1 

43.2 

K s band 

6(R.A.) 

42.5 10.8 

7.1 

4.1 

44.6 

6(Dec.) 

41 11.0 

7.1 

4.1 

43.2 


R.A. (Bulge) 


Dec. (Bulge) 


11:32:34.8515 ±0.0051 

+53:04:04.475 + 0.044 


R.A. (Red blob) 


Dec. (Red blob) 


11:32:34.8550 ±0.0051 

+53:04:04.512 + 0.044 


R.A. (Radio) 


Dec. (Radio) 


11:32:34.8497 ±0.0005* 

+53:04:04.527 + 0.004 

Notes. Upper table: Total astrometric uncertainties involved 

in our re- 


suits interpretation for each band. Columns: 1) Total intrinsic SDSS 
astrometric uncertainties, 2) z-J bands matching uncertainties during 
the SUBARU astrometric calibration, 3) centering and/or shifting un¬ 
certainties during J,H,K S centering and the J-H,J-K S alignments, 4) 
intrinsic radio map astrometric uncertainties, 5) quadratic addition of 
all the above. Lower table: From top to bottom, positions of the peak 
flux of a) the stellar bulge of NGC3718, b) the offset red blob and c) the 
offset e-Merlin 6 cm radio emission. 

w This is the position of the peak of the radio emis sion of NGC3718, 
using the position of J1 146+5356 quoted in section |22| Assuming the 
position used in earlier e-MERLIN observations, our R.A. and Dec. 
measurements namely, 11:32:34.8534 + 0.0005 and +53:04:04.523 + 
0.004, are, within the uncertainties, in very good agreement with the 
position published by |Krips et al.| ( |2007] ). 


4. NIR color maps 

Having aligned and calibrated our NIR images, we proceed by 
producing the J-H and H-Ks NIR color maps. NIR color maps 
are an important tool as they provide vital information about the 
color distribution of a galaxy, such as whether reddening affects 
the colors and whether light comes purely from an “ordinary” 
stellar component or if it contains a contribution from a nuclear 

9 We magnify a 51 x 51 pixels 2 (FoV 20 x 20 arcsec 2 ) part of the 
SDSS z band centered on the nucleus of NGC3718, by a factor of ~ 40, 
equal to the scale ratio of the SDSS and SUBARU data (namely ~ 0.4 
and ~ 0.01 respectively). The new, 2040 x 2040 pixel 2 , SDSS z 
band replica is then aligned with our SUBARU J(H and K s , trimmed to 
2040 x 2040 pixel 2 ) band(s), following the method described in section 
3J"1 
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component as well ( [Glass & Moorw ood 1985 ). The colo rs of 
“ordi nary” galaxies are, J - H = 0.78 and H - K = 0.22 ( [Glass 


|1984| ). Light captured by the J band is produced mainly from 
stars, whereas in K s , light can come from both stars and hot 
dust. A black body with a temperature below the sublimation 
temperature of dust (~ 1500 K), will radiate at A > 1.9/tm, 
mainly affecting the K s band but, potentially, also the H band 
(Acentral = 1.63/im) to some extent. In this context, J-H is 
mostly indicative of extinction, as dust is already evaporated, 
and is acting as a scattering source, while H - Ks, is indicative 
for hot dust emission. The NIR color maps are shown in figure 


H-Ks 



0.64 

0.60 

0.56 

0.52 

0.48 

0.44 

0.40 

0.36 

0.32 



02.44 

35.080 34.861 34.641 


R.A. (ll:32:ss.sss) 

Figure 8. H - K s (top), J - //(bottom) NIR color maps. 


Our J-H map is fairly symmetric. An approximately el¬ 
liptical red excess can be seen surrounding the nuclear region 
extending from SE to NW. This could be interpreted as local 
extinction, since there is no apparent large scale constant color 
gradient , which would indicate that the dust lane (lying SW and 
outside the borders of the light distribution, see Appendix[A]) ex¬ 
tends all the way to the projection of the center of the bulge. The 
H - Ks map also shows a similar, but better defined, elliptical 


red excess around the center (sharing the same orientation with 
the excess in the / - H map), indicative for the presence of hot 
dust emission. 

The overall symmetry of the NIR color maps can be more 
clearly seen on the radial profile plots of figure [9] In the four 
different directions of measurement (see caption of figure [9]), the 
plots do not show any sign of a large scale constant color gra¬ 
dient, which would appear as an asymmetry with respect to the 
center and it would indicate the presence of foreground dus{^ 
The only asymmetry that can be seen is between the blue and rec 
profiles outside of ~ 0.5 arcsec. This is due to the fact that the 
directions of measurement are perpendicular to each other and 
cross the center along the “major” and the “minor” axis of the 
red elliptical excess respectively. In any case, both of the com¬ 
ponents are highly symmetrical with respect to the center. 


H-Ks vs Re I. offset 




-1.0 -0.5 0.0 0.5 1.0 

Relative offset (in arcsec) 


J-H vs Rel. offset 



-1.0 -0.5 0.0 0.5 1.0 



Relative offset (in arcsec) 


Figure 9. Radial profiles of the H - K s (left column) and J-H 
(right column) NIR color maps, derived from two sets of perpen¬ 
dicular directions measured counterclockwise from the North. 
These are, 120°/210° (blue/red points respectively) in the up¬ 
per row and 170°/260° (green/magenta points respectively) 
in the lower row. Negative—positive relative offset reflects a 
S(E,W)->N(W,E) direction. 


The most interesting feature of the H-Ks map, however, is 
that the center of the red blob (i.e. the peak of the H - K s map) 
lies slightly offset from the center of the stellar bulge, by ~ 50 
mas. For this to be more clearly seen, we produce the symmetri¬ 
cal H - K s map (using the H and K s images of figure [7]) and we 
over-plot its peak contours (i.e. the center of the stellar bulge), 
alongside the peak contours of our normal H-Ks map, in fig¬ 
ure [TO] We also manage to extract the same offset length, along 
the same orientation, by applying a high pass filter directly on 


10 A non-uniform, progressively thinner, dust extension northward of 
the dust lane, would translate into progressively excessive dust towards 
the S/SW in comparison to the N/NE. This would lead to excessive red¬ 
dening, resulting in an asymmetry in the figure [9] plots, in the form of: 
progressively more dust —> progressively increasing reddening —> pro¬ 
gressively fainter /(//) —> progressively larger J - H (H - K s ) towards 
the S/SW. A uniformly distributed foreground dust extension, however, 
could be present, but in this case the reddening should be uniform as 
well, a fact that would not affect (at least) the (qualitative) interpreta¬ 
tion of our results. 
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the reduced K s image, a result “free” from any geometric trans¬ 
formation uncertainties. The corresponding central astronomical 
positions for the stellar bulge and the offset red blob, respec¬ 
tively, are shown in table [2] 


H-Ks _ 

-| 0.64 


0.60 



R.A. (11:32:34.sss) 
Ks band 



R.A. (11:32:34.sss) 

Figure 10. Top: H - Ks closeup. The red blob contours (red) are 
offset from the symmetrical light image’s peak contours (i.e. the 
center of the stellar bulge, in blue). Bottom: K s band high pass 
filtered image closeup. The peak flux is offset from the photocen¬ 
ter of the K s band. In both images, the red circle indicates the 
position of the peak value of the offset red blob, while the blue 
circle indicates the position of the center of the stellar bulge, 
as derived from our contour analysis. Their radii indicate the 
positional uncertainty of each position, which is ~ 7 mas. The 
yellow(top)/cyan(bottom) line is the length of the offset, namely 
~ 50 mas. 


4.1. Light decomposition 


Following Glass & Moorwood ( j!985| , we use their figure 7 
(see caption of figure ITT] for an explanation), of theoretical 


loci for mixtures of “ordinary” galaxy colors and blackbodies 
of various temperatures constrained on the color-color diagram 
of their sample, in order to decompose the light of NGC3718. 
An azimuthal average, binned in 5 equally spaced (~ 50 mas 
step) radius bins, of the innermost 0.5 arcsec of the nucleus of 
NGC3718, can be seen in (the zoomed-in section of) figure [TT] 
A decline, from a 50% stellar - 50% dust (at Tioo ~ 1000 K) 
light mixture in the innermost ~ 100 mas, to a 60% stellar - 40% 
dust (at T 2 oo ~ 800 K) in the innermost ~ 200 mas, which drops 
even further to a 65% stellar - 35% dust (at 7300-500 ~ 700 - 500 
K) in the innermost ~ 300 - 500 mas, can be clearly seen in 
this plot. This indicates the presence of significant contribution 
from hot dust in the central region of NGC3718, suggesting an 
environment ideal for SMBH accretion. 



0 1 2 
H - K 


Figure 11. 


Glass & Moorwood 


(1985) two color diagram. 


Moving away from the “ordinary” galaxy colors (black circle 
at J - H = 0.78 and H - K = 0.22) and along the horizontal 
lines in the loci of mixtures, there is an increasing percentage 
(10% per vertical line) of contribution from hot dust emission, 
while moving along the vertical lines, the temperature of this 
dust component increases (T ~ 500 - 1500 K). 


4.2. Over-plotting the e-Merlin radio map 

The high angular resolution (~ 100 mas) 6 cm radio e-Merlin 
map, with a pointing accuracy of ~ 4 mas (see section |T2| ), is 
shown in figure [12] The outer contours are set to ~ 4<x sigma 
above the background, so we consider the flux density enclosed 
within these contours to be a true detection of extended radio 
emission, at least, in the position angle shown by the red dashed 
line in figure [l2| 

Over-plotting the radio map on our H - Ks NIR map, re¬ 
veals that, the radio emission appears to originate from a posi¬ 
tion closer to the offset red blob than to the center of the stellar 
bulge (figure [L3]). This is yet another very interesting result, be¬ 
cause it could indicate that the SMBH does not lie exactly at 
the center of the bulge, but it is offset by ~ 50 mas. This, at a 
distance of ~ 17.4 Mpc, yields a projected offset of ~ 4.25 pc. 
The presence of radio emission is expected, since NGC3718 is 
a LLAGN, so the closer proximity of the radio emission to the 
offset red blob, implies that the SMBH accretes matter normally, 
but it does so slightly away from the center of the stellar bulge. 
The coordinate uncertainties however (see table [2]), are of the 
order of the offset, though little smaller. This suggests, that a po¬ 
tential coincidence of the radio emission with the offset red blob, 
is reasonably likely to be true. The astronomical position of the 
peak radio flux, is shown in table [2] 
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Merlin 6 cm 


0.0056 

0.0048 

0.0040 

0.0032 

0.0024 


0.0016 

0.0008 


0.0000 

'“.905 .851 .797 

R.A. (11:32:34.sss) 

Figure 12. e-Merlin 6 cm radio map. The outer contours are set 
to ~ 4<x above the background where the flux is ~ 2 orders of 
magnitude less than the peak flux. The third inner contour is 
~ 40 - 50% of the peak flux. The red and blue dashed lines, 
represent the position angle of the SE-NW and the NE-SW can¬ 
didate small scale bipolar extended structures, respectively. The 
beam is shown in the upper left comer. 



H-Ks 



0.64 

0.60 
| 0.56 
0.52 

0.48 

0.44 

0.40 

0.36 

0.32 


1 


Figure 13. Same as the upper image of figure 10 but with the 
radio contours over-plotted (in the position quoted in table 
The same image with the equivalent to the Krips et al. (2007 
position, is shown in Appendix [C] figure [Cl] ). The blue and red 
circle radii, indicate the positional uncertainty of each position, 
namely ~ 7 mas. The black dashed ellipse indicates the total 
coordinates uncertainty of table [2] 


Spatial offsets, potentially coincident with AGN offsets, are 
predicted by current theory. A possible interpretation could be 
that this is the ca se of a SMBH recoil (e.g. Komossa 2012| 
Sija cki et al.|2011| ). In short, in the case of a merger, the SMBHs 
present in the centers of the merging galaxies will form a bi¬ 
nary ( [Begelman et al.|1980| ) which will eventually coalesce due 


to anisotropic emission of gravitational waves (e.g. Bekenstein 
[1973| ). We need, however, more evidence to favor such a view of 
NGC3718, and will resume the discussion on this subject later 
on this paper. 

We also confirm the presence of the ~ 0.5 arcsec NW can¬ 
didate jet component, that |Krips et al.| ( [2007] ) see in their Merlin 
18 cm map. Moreover, we detect a diametrically opposite placed 
~0.5 arcsec extension towards SE. This could be considered as 
a possible counterpart of that jet, since, both lie along the same 
orientation and are of the same length. Finally, a less extended 
but brighter structure towards SVyj’’] and, a fainter tail-like one 
towards NE, both of - 0.3 arcsec in length and sharing the same 
orientation, could be, potentially, indicative of a small scale X- 
shape radio structure, often associated with “spin-flip” processes 
in SMBH recoils (e.g. Liu et al.|2012| ) or with the two jet emit¬ 
ting members of a close SMBH binary ( |Lal & Rao 2007]). 

Regarding the 6 cm e-Merlin radio flux of NGC3718, we 
measure a peak flux of 8.85 ± 0.07 mJy and an integrated flux of 
9.97 ± 0.08 mJy, both of w hich are larger than the values pub¬ 
lished by [Krips et al.| ( [2007] ), namely, 5.3 ±0.1 mJy and 6.1 ±0.3 
mJy, respectively. This reflects an increase in peak flux of - 67% 
and in integrated flux of - 63%, indicating that NGC3718 is 
a variable radio source, a picture consistent with an accreting 
SMBH. 


5. Scaling relations & classification 

5.1. To "bulge" or to "pseudo-bulge”? 

The most famous scaling relations between a host galaxy and its 
SMBH, are the Mbh - L and Mbh - cr , relating the luminosity 
and the stellar velocity dispersion of the bulge, respectively, with 
the mass of the SMBH. 

We can estimate the absolute magnitude M Ks ^ u i ge of 
NGC3718 from our data. Given the fact that we do not per¬ 
form any decomposition into different components (disk, bulge 
etc.), the (contaminated) bulge magnitude we derive, should pro¬ 
vide an upper limit for the estimation of the Mbh- We perform 
aperture photometry on our K s band images (figures [3] and [7} 
of NGC3718, using three different aperture sizes, from 3 - 5<x 
above the background, with an average aperture of r apeK - 3 
arcsec]^] We average the flux of those measurements and we use 
a (rather large) 0.2 mag assumecp] error. The K s band absolute 
magnitude of NGC3718 is: 

Mk s ,bulge = -21.1 ± 0.2 mag . 

For the velocity dispersion however, we will need spectro¬ 
scopic information. Luckily, NGC3718 is a very well studied ob¬ 
ject. Ho et al. ( |1997a| ), in their Palomar AGN spectroscopic sur¬ 
vey, were able to extract several useful parameters, among those, 
the FWHM[NII] = 371 km s -1 line. Assuming a Gaussian dis¬ 
tribution and using FWHM = 2\/(2ln(2))cr « 2.355<x, we can 
get an estimation for the velocity dispersion, of cr « 157 km 
s -1 . However, [W//]/U6583A is considered to be a better tracer 
for the ionized gas motion. A number of authors have suggested 


11 In the position angle of the blue dashed line, which could be identi¬ 
fied as the westward candidate jet component, [Krips et al.[ ( |2007| see in 
their EVN 6 cm map. The same structure, aligns also with the jet-like 
structure seen at arcmin scales by the VLA at 1.49GHz in C/D-array 
configuration by Condon|(|1987|). 

12 The average difference of the flux measurements on the images of 
figure [3] from those of figure^] is ~ 0.1 mag which we consider negli¬ 
gible. 

13 To account for the statistical deviation of the measurements, for 
systematic uncertainties in the calibration process etc. 
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that the gas may rotate at a different speed t h an the star s (e.g. 
Pig natelli et al.| (200 TT ), |Caldwell et al.| ( |l986| )). |Ho| ( |2009| ) com¬ 
bined gas and stellar velocity dispersion measurements an d they 
conclude that <x* and cr g correlate well. Ho et al.| (J2009J) calcu¬ 
late a velocity dispersion for NGC3718 of cr = 158.1 ± 9.6 km 
s"\ using the Ca + Fe absorption feature at 6495A, value which 
we adopt. This suggests a ^ ~ 1 for NGC3718. The author 
argues that “...as the gas derives principally from mass loss from 
bulge stars, its kinematics should generally track the kinematics 
of the stars. But because the gas is collisional and experiences 
hydro dynamical drag against the surrounding hot medium, we 
expect it to be kinematically slightly colder than the stars. In the 
absence of additional energy input from other sources, we an¬ 
ticipate < 1 , as observed. As additional energy is injected 
into the system, for example from activation of the central black 
hole, the gas gains energy, to the point that cr gas approaches or 
even overtakes crf\ This picture is consistent with the fact that 
NGC3718 is indeed an LLAGN. This could be considered as 
an indirect indication that, NGC3718’s central engine provides 
enough AGN feedback, to heat up the gas and cause this strong 
agreement between cr gas and <r*, information which will be used 
later on. 

Having an estimation of both of these quantities, we use the 
M B h - M Ks ,bulge an d 4 ^bh - cr relations from|Kormendy & Ho 
( |2013| ), namely: 


projection effect of the prominent dust lane of NGC3718, which 
is the result of looking through a disk/ring of dusty gas on po¬ 
lar orbit around the, nearly face-on, stellar disk. They also report 
that the latter must be substantially free of cool gas, since they 
do not detect any HI emission on the plane of rotation of the stel¬ 
lar disc, a characteristic of lenticular galaxies. Additionally, the 
LI.9 spectral classification (Ho et al. |1997b| , means that broad 
line H a emission is detected, suggesting that our line of sight 
provides a direct view towards the nucleus of the host galaxy 


which, consequently, has to be closer to face-on. [Sparke et al. 
( |2009[ ) conclude that NGC3718 is a Polar Ring Galaxy, a rare 
class of objects, often hosting a lenticular galaxy, which as an 
early-type object, shares many kinematic and other properties 
with elliptical galaxies. Our data also favor this picture, since 
we do not see any traces of a nuclear bar. The elliptical region 
around the nucleus (H - Ks map of figure [8]), looks more like an 
almost face-on disk rather than a bar. 

Moreover, NGC3718 does not fulfill at least three of the cri¬ 
teria of |Kormendy & Kennicutt (2004), for classifying a bulge 
as a pseudo-bulge, namely: 

1. Following |Kormendy & Bender| ( |2013| ), we calculate the 
Faber-Jackson correlation for coreless and core ellipticals re¬ 
spectively, using Ly = 1.20 x 10 10 L o , namely: 


Ly 


log 


4^bh 

10 9 M o 

4^bh 


= -(0.265±0.050)-(0.488±0.033)(M^^+24.21)(5) 


10 n L Q 


10 n L o 


(0.67+0, 


= (0.79+0.11) 


^ ^\25tikms 1 ) 

( * \ 

\250kms~ l ) 


3.74+0.21 


8 . 33 + 1.24 


I58kms~ l (7) 


199 kms~ l (8) 


108 = -» 5 09 ± 0 fl 49) + <4.384 ± 0.287)log( 5s £ FT )« i , 

Equation [5] suggests an SMBH mass: 

M MsH-M KsMlge = j 65+1.64 x 1Q 7 Mo ^ 

whereas equation^] suggests a mass: 

A& _<r = 1.11+1'“ x 10 8 Mq _ 

Equations [5] and[6| despite the fact that they are derived from 
the same sample by Kormendy & Ho ([2013), differ by almost 
an order of magnitude in terms of M BH . Adopting M^ H ~ (T , as a 
more robust result and calculating the magnitude that the bulge 
should have to justify this mass, we find a minimum difference 
of: 

SM Ks -1.7 mag 

or, in other words, NGC3718 should be at least - 4.8 times 
brighter. 

But what can cause such a disagreement between the two 
scaling relations, in the case of NGC3718? A first approach 
could be that, this is not a bulge but a pseudo-bulg^J so these 
equations do not apply ( [Kormendy et al.||2011| ). NGC3718 is 
currently classified a s an SB(s)a pec, m eaning that it has a bar 
and, according to [Kormendy & Kennicutt ( 2004| , barred galax¬ 
ies preferably have pseudo-bulges. 

On the contrary, several au thors (e.g |Pott et al.|20 04; Krips 
et al.||2005{ |Sparke et al.||2009| ) have studied the gas dynam- 


NGC3718’s kinematic properties (Jog & Chitre 200Z| im¬ 
ply the presence of a random motion supported bulge and 
its velocity dispersion (Ho et al. j2009 ) seems to be in good 
agreement with both versions of the Faber-Jackson correla¬ 
tion (given their large scatter), although a closer match to 
that of a coreless elliptical (equation [7]). In any case, the ob¬ 
served cr is not much smaller than the cr predicted by these 
relations, which would be indicative of a rotationally sup¬ 
ported pseudo-bulge. As pointed out by |Kormendy & Ho 
(2013), “Classical bulges are essentially equivalent to core¬ 


less ellipticals ”, a picture consistent with the classification of 
NGC3718 as a lenticular galaxy. Considering also, 

2 . that we do not detect any signs of a nuclear bar and 

3. that there is no sign of significant star formatior|^ ( Sparke 
[et al.|2009[ ), 

allowing us to assume that NGC3718 behaves as an elliptical 
galaxy and, therefore, it should follow the host galaxy - SMBH 
mass scaling relations. But why do the M B h - M Ks ,bulge an d 
M B h - cr relations point towards different M B h ? 


ics and have successfully fitted tilted rings on NGC3718. They 
all agree that the gas orbits are nearly edge-on near the center. 
|Sparke et al.| |2009| suggest that the classification as a barred 
galaxy is misleading. They see the apparent nuclear bar, as a 


5.2. Light “deficit” or cr “surplus”? 

Given the fact that equations [5] and |6| are calibra ted against the 
same galaxy sample by |Kormendy & Ho ( 2013| , we can safely 
assume that t hey should, more or less, agree. As previously 
noted (section HT we consider cr to be more accurately de¬ 
termined than M Ks £ U i ge , since, as it is measured from two dif¬ 
ferent linewidths and the Faber-Jackson relation, shows surpris¬ 
ingly good consistency. 


14 A pseudo-bulge is a central high density stellar region similar to a 
classical bulge, but with properties closer to disk-like objects (e.g. sup¬ 
ported mainly from rotation rather than from pressure due to random 
motions of stars). Their formation mechanisms are believed to be dif¬ 
ferent than the ones that formed classical bulges. 


15 We support the view of Sparke et al. (2009 ). We examine GALEX 
near and far UY images and we find that, NGC3718’s UV flux is - 1 
order of magnitude less than the flux of the starburst galaxy NGC4449, 
which indicates some amount of star formation, but at the same time, 
this is not high enough to be characterized as “starburst”. 
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Our first choice therefore, is to examine the robustness of 
Mk s , bulge- In order to better understand our light distribution, 
we produce the surface brightness profile of NGC3718, from a 
joined data set consisting of our SUBARTp|and 2MASSp^l K s 
data, against the radius. For this, we use IRAF’s task ellipse once 
again, which returns an azimuthally averaged surface brig htness 
profile, as shown in figure [14] Many authors (e.g. Kormendy 
|et al.|2009| argue that, fitting a Sersic profile to larger radii of 
an elliptical galaxy and then extrapolating it inwardly, is an ac¬ 
curate way to describe the surface brightness profile. 

Following this scheme, we fit our combined SUBARU and 
2MASS K s data set of NGC3718. Our best fit consists of a Sersic 
profile for the bulge with an index n = 4.85, combined with a 
dual exponential disk component with scale lengths R dlskl ~ 7 
arcsec and R dlsk2 ~ 60 arcsec respectively, as presented in figure 


14 We choose the lower limit of our fitting range (~ 50 arcsec), 
based on the uncertainties of the 2MASS data points (at ~ 50 
arcsec the flux is ~ 2 cr above the background). The upper limit 
is then determined through ^ 2 minimization for that given lower 
limit, in the range between the (candidate) upper limit and ~ 20 
arcsec, in order to minimize the fit with respect to the higher 
signal-to-noise data points (at ~ 20 arcsec the flux is ~ 6 cr above 
the background). Our best fit has a x 2 ~ 2, for a data range 
between ~ 1.5-50 arcsec (~ 500 data points). Our modeled 

bulge magnitude estimation is Mk s ,bulge -22.85 magnitudes 

(m Ks ,bulge ~ 8.35). When put into equation [5] this “recovered” 
magnitude, suggests an Mbh of the order of: 

M bh KsMlge = 1-18^59 X 10 8 M o , 

which, surprisingly enough, covers the gap between the 
M B h ~ L and M B h - cr relations (section [571] ), almost perfectly. 

A similar result, was previously noted by |Dong & De| 
Robertis[ ( |2006[ ), who used 2MASS data of NGC3718 along with 


the bulge/disk decomposition algorithm GALFIT, in order to 
separate the bulge from the disk and obtain magnitudes for both. 
Their best fit consists of a Sersic bulge profile of index n = 5.5 
with M K°7uige ~ - 22 ; 75 magnitudes (m°°J ulge ~ 8.45) and a disk 
with scale length Rf sk ~ 60 arcsec. This magnitude estimation 
is almost identical to ours. In this case, equation [5] then gives: 

m2 H ~ KT ^ = 1.05!^ x 10 8 M o . 

What both of these fits reveal is very interesting, since, in¬ 
stead of a Sersic profile near the center, we get a plateau. The 
fit, describes the surface brightness profile very well at radii 
> 1-1.5 arcsec, but it “fails” near the center, where the model 
suggests that an amount of “missing light” is present, seenas the 
difference between the fitted and the observed curves. Kormendy 
et al. ( 2009| suggest that a surface brightness profile of this kind, 


is a characteristic seen in core galaxies. 

In short, a “core” is an elliptical galaxy whose measured sur¬ 
face brightness profile shows a break towards a shallower light 
profile when compared to an inward extrapolation of its outer 
Sersic fit, revealing an amount of light “deficit”. The opposite 
stands for coreless galaxies. These are objects whose measured 
surface brightness profiles show a break towards a steeper light 
profile with respect to their, inwardly extrapolated, outer Sersic 
fit, revealing an amount of light “surplus”. 

So, as our analysis suggests, the light “deficit” is, most 
likely, the source of the disagreement between the M B h - L and 


16 Sample in the inner ~ 2.5 arcsec. 

17 Sample between ~ 2.5 to ~ 70 arcsec. 



Figure 14. First figure: Surface brightness in the Ks band ver¬ 
sus log(R) in pc. The blue dashed line represents the inwardly 
extrapolated Sersic profile of the bulge, while the black dashed 
line represents the dual exponential disk profile. The blue points 
represent the 2 MASS Ks data, while red points represent our 
SUBARU K s data. Second figure: Same as the first figure, but 
with the different components added. Here the red dashed line 
represents the fitted model whereas the black points represent the 
complete joined data set. Third figure: Surface brightness in the 
K s band versus R (in arcsec) of the innermost ~ 3 arcsec. The 
red dashed line represents the fitted model, while black points 
represent (up to ~ 2.5 arcsec) our SUBARU K s data. Fourth fig¬ 
ure: Same as the third figure, but plotted for the complete data 
range (~ 70 arcsec). The green lines represent the uncertainties 
of the data points. 


M B h - cr relations 18 19 As discussed by several authors (e.g 


18 We have to note at this point, that we agree with Hopkins et al. 
42008) on the “danger” of doing three component fits. Our fit is indeed 
very sensitive to the adopted lower and upper limits. This is the reason 
why we resort in trusting the statistical robustness of our sample and set 
statistical, yet as physical a s possible, constrains. The fa ct that we con¬ 
firm, more or less, the fit of Dong & De Robertis ( 2006) indicates some 
consistency. Moreover, we know that light is “missing” from the central 
region, since the published aperture photometry values for NGC3718 on 
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Kormen dy et al. | |20Q9t |Kormendy & Ho 2013; Gual andris & 


Merritt 2008), such a light “deficit’ | zu | could be indicative of co¬ 


evolution driven by a combination of SMBH binaries followed 
by SMBH recoils. The light “deficit”, can be retrieved by the dif¬ 
ference between the observed and the expected brightness (fig- 

2.2 x 10 9 Mr, 


ure 


14). In our case, it translates to an M^ C3718 


and a ~ 19. When put onto the M^ e f versus M B h plot of 
|Kormendy & Ho| ( j2013| (their figure 30), it falls on the left of 
and away from the main distribution (M^f ~ 5M BH ), but well 
within the observed upper limit of M^ef ~ 50 M BH . Moreover, 
if NGC3718 is viewed as an advanced merger ( (Chitre & Jog[ 
|2002[ ), then it falls much closer to the main distribution when 
compared to NGC4486B, the only other ongoing merger on this 

m NGC44S6B 

P lot > with a 

a mass “deficit” of < 


(e.g. 
from 

MNGCtttk 

def 


Hopkins et al 
Schwarz (198 


l »: 


> 50. Additionally, for a real(istic) core, 
1 % of the total galaxy mass is expected 


2008). Using M\ 


NGC3718 _ 400 X 10 9 M o 


and calculating the percentage, we get an 
0.5 %M^p C3718 , which is well within the realistic 

ap- 


limits. These facts, make NGC3718’s relatively large M^f 
pear as, at least, physically plausible. 

Finally, our fit is consistent with previous studies. Chitre 
& Jog| ( |2002| and |Jog & Chitre) ( |2002| studied 2MASS and 


kinematic data of NGC3718, for which they find elliptical-like 
kinematics mixed with spiral-like photometry. Bournaud et al. 


( |2005| ), attributes the “transitional” nature of these objects to re¬ 
sult from mergers with mass ratios > (3 - 4.5) : 1, remnants 
which closely resemble the properties of lenticular galaxies. The 
fact that we and [Dong & De Robertis ( 2006), see both the expo¬ 
nential disk and the random motion supported bulge in the sur¬ 
face brightness profile means that NGC3718 should be classified 
as a lenticular galaxy which, consequently, strongly implies that 
NGC3718 is a merger remnant. 

The mass “deficit” and the indications for a merger, are com¬ 
patible with the evidence we have so far, that NGC3718 might 
host an actual SMBH recoil. We discuss further this subject in 
the discussion (section [7]). For now, we adopt the average value 
of the equations [5] (ours and Dong & De Robertis] ( |2006[ )) and [6] 
for the rest of our calculations, as a good estimation for the M BH , 
namely: 

M M?C3718 = 1_ 11 + 1.09 x 10 8 M _ 


± BH 


0.55 


5.3. Core or core less? 

As we have seen up to now, NGC3718 appears to show proper¬ 
ties that belong to both core and coreless ellipticals. For exam- 


NED are closer to our SUBARU measurement than to the “recovered” 
magnitudes. We can not, however, rule out the case of a bit different 
Sersic index with a small “extra light” component near the position of 
the transition from the disk to the bulge dominated component (log(R) 
~ 2 - 2.5 pc in figure plj). Such a transition would be very difficult, if 
not impossible, to be distinguished in this case due to the presence of 
the disk. Provided that NGC3718 is a merger remnant, the low star for¬ 
mation activity is expected to have created an amount of post merging 
stars, that could, in principal, have resembled an “extra light” compo¬ 
nent. We think that this would not change dramatically our overall view, 
especially regarding the amount of the “missing light” which seems to 
be quite dominant. 

19 The sensitivit y issu e on the adopted upper and lower limits, is also 
mentioned by |Jog & Chitre| ( [2002| ) and they decide not to make such a 
fit. Our SUBARU data, however, resolve the innermost ~ 2.5 arcsec of 
NGC3718 much better, providing us with a richer data sample to work 
with. 

20 Observed up to now mostly on bright giant ellipticals. 


pie, its velocity dispersion fits more to that of a coreless galaxy, 
whereas the presence of “missing light” is a characteristic usu¬ 
ally seen in core galaxies. It is reasonable, therefore, to inves¬ 
tigate further similarities and differences with either of these 
categories. The classification criteria of the E-E dichotom)p" 
have been discussed by many authors (e.gjEmsellem et al.|2007 


Tho mas et al.|2005 Cappellari et al. 2007J) and the most updatec 
classification criteria for core and coreless elliptical galaxies, are 
presented by |Kormendy | ( [2009] ). 

Supporting a classification as a coreless object is the visual 
magnitude My GC371s = -20.73 mag, which is fainter by almost 

an order of magnitude than the My -21.5 mag, the “border’ p 7 ] 

of the E-E dichotomy, lying on the bright end of the coreless side 
(see figure 2 of [Kormendy |2009| ). Additionally, the presence of 
weak/not dominant radio emission, as we show in the e-Merlin 
6 cm map (figure [T2|) and as measured by Krips et al. ( |2007| and 
by us (see sectio n|4.2|), which, along with NGC3718’s kinematic 
properties ( [Jog & Chitre|2002] ), are characteristics usually seen 
in coreless objects. 



Figure 15. A4 parameter versus radius for NGC3718. From top 
to bottom, A4 parameter for the K s (red/top), H (blue/middle) 
and J (green/bottom) bands respectively, versus R in arcsec. In 
all plots, the dashed black line represents the border between a 
“disky” and a “boxy” object. 


On the other hand, the most obvious observational evidence 
to support a classification as a core object, is the “missing light” 
that we observe in the surface brightness profile. The fact that 
the Sersic profile has an index of n > 4 and the presence of 
moderate to low hard X-ray emission (L^ 7B y T s ~ BAT = 6.46x 10 41 

ergs sec -1 from Tueller et al. (2010)), may further strengthen a 


classification as a core object. 

Regarding the isophote contours shape of the images of fig¬ 
ure]^ we see the same general trend in all bands, as we present in 


21 Elliptical galaxies form two distinct categories, each one defined by 
specific fundamental properties such as core or coreless, boxy or disky 
etc. 

22 The E-E dichotomy “border” is calculated for Hq = 72 km s -1 . For 
this value the absolute magnitude of NGC3718 becomes My GC31ls = 
-20.60 mag. 
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figure [l5] The A4 parametef^is mostly negative up to a radius 
of ~ 1-1.5 arcsec (essentially, the radius of the core we ob¬ 
serve), consistent with a boxy (core) object, while further out it 
steadily raises above zero and drops again, to remain on the pos¬ 
itive side, up to a radius of ~ 3 arcsec, consistent with a disky 
(coreless) object. This could be indicative of an inside-out trans¬ 
formation of NGC3718 from a coreless to a core galaxy. The 
available classification criteria, are summarized in table [3] 


Table 3. E-E dichotomy parameters. 


Criteria 

core 

coreless 

Reference 

Visual Magnitude 


/ 

1 

Surface Brightness profile 

/ 


5 

Isophote Contours shape 

/ 

/ 

4 

Sersic index 

/ 


5,6 

Radio emission 


/ 

3,7,8,9 

X-ray emission 

/ 


2 

Kinematics 


/ 

10 


Notes. E-E dichotomy classification parameters for NGC3718. 


References. (1) NED; (2)|Tueller et al.|(|2010j); (3) Figure[l2l (4) Figure 
fl5l (5) Figurepl (6)|Dong & De Robertis C2O05J; (7)|Pott et al.|p004); 
(8) |Krips et al.| ( |2005| K(9) lKrips et al.| ( |2007| >; (10) |Jog & Chitre||2002] );. 


The picture we get from the above discussion is that 
NGC3718 is most likely an intermediate (between core and core¬ 
less) galaxy, a fact that makes the attempt for a possible inter¬ 
pretation of the results more complex. It provides, however, an 
additional piece of evidence in favor of a view of NGC3718, as 
the “smoking gun” of a past/ongoing merger, in the sense that 
the E-E dichotomy refers to early-type objects, believed to have 
formed in mergers. 


6. Electromagnetic signatures of an SMBH recoil 


So far, the observations imply that NGC3718 is a merger rem¬ 
nant, with a core in its surface brightness profile, combined with 
a spatial and a, potentially coincident, AGN offset, all of which 
are compatible with SMBH recoils. 


If this is indeed the case, according to |Komossa| ( 2012| , we 
should expect to see some characteristic electromagnetic signa¬ 
tures. The most prominent of those, would be a velocity shift 
between the broad and the narrow component of a given line. 
Briefly, the idea is that a region of the order of size of the Broad 
Line Region (BLR) should be gravitationally bound to the re¬ 
coiled SMBH, resulting in a relative motion of the BLR with re¬ 


spect to the, much larger, Narrow Line Region (NLR). Ho et al 
( |1997b| ) fit a narrow and a broad component to the Ha line and 
they find that the centroid of the broad Ha is redshifted with re¬ 
spect to its narrow counterpart by 430 km s -1 (their figure lOe). 
If we also account for systematic offsets of the order of 50 - 100 
km s“\ between broad and narrow line components (Peterson 


1997| ), we adopt a kick velocity of the order of: 

Vkick ~ 355 km s -1 


23 A4 Fourier coefficient indicating the deviation from an ellipse. 
“Boxy” objects have a A4 < 0, whereas “disky” objects have a A4 > 0. 


Based on the Vkick and M ^ C3718 estimations, we can calcu¬ 
late the gravitationally bound radius, according to the relation 
( [Merritt et al.|2006[|Komossa|2012| ): 


Rbound 


/ G x M bh 

\ Vkick ^ 


0.4 x 


Mi 


BH 


\ 10 8 M C 


\ / Vkick \ ^ 

) / \ 10 3 kms~ l ) P C 


(9) 


which results to an Rbound ~ 3.53 pc. Since the gravitationally 
bound radius is inversely proportional to the itoTd and due to 
the large uncertainties of the M^ C3718 , this should be viewed 
as an indicative, yet relatively uncertain, estimation. With this in 


mind, we can estimate the size of the BLR, using the relation 
from |Zhang et aL| ( |2007| ) : 


Rblr = qx 


G x Mbh 


FWHM Ha , Broad 1 


( 10 ) 


which, depending on the adopted value of q, gives an R B lr ~ 
0.06 - 0.12 pc (q = 0.75 - 1.33), usin g FWHM Ha , Broad * 2350 
km s " 1 for NGC3718 from Ho et al. (1997a). One last quantity 


to estimate, is the sphere of influence of the SMBH, given by the 
relation: 


Ri 


infl. 


_ | 6 x M bh j 


(ID 


which equals to R m fi. ~ 19.1 pc for the values we have adopted 
up to now in this paper. The latter is ~ 5 times smaller than the 
core radius, i.e. ~ 100 pc. 

These values are useful, in order to better understand the 
nuclear dynamics of the system. Specifically, the R blr and the 
Rbound > strongly suggest that a BLR of the predicted size would 
likely be, entirely, gravitationally bound to a recoiled SMBH. 

Lor our final calculation, we use the images of figure [7] and 
we measure values below the resolution limit. The fact that these 
images have not been extensively tested for their scientific ac¬ 
curacy, makes them useful only for rough (order of magnitude) 
estimations (see section [33] ). 

With this in mind, the presence of a Hyper Compact Stellar 
System (HCSS) is expected around an SMBH recoil. According 
to |Komossa| ( |2012| ), an HCSS should be found within an Rbound 
and its total luminosity should be of the order of the luminosity 
of a globular cluster. We attempt to measure this, by measuring 
the Ks band fluxes of figures [7] and [3] at the position of the off¬ 
set red blob, within an Rbound and taking the difference. These 
two images differ only by the presence of the offset red blob in 
that region. The rest of the emission through this column, should 
be a good approximation of the “foreground” and “background” 
emitting sources, respectively, with respect to the center. This 
luminosity difference is: 

Lhcss ~ 3.5 x 10 4 L o , 

whereas this converted to mass (see section [I]), yields a: 

Mhcss ~ 2.7 x 10 4 M o , 

which is within the predicted luminosity and mass range for a 
globular cluster. 

These findings, further strengthen the possibility that an 
SMBH recoil is indeed present in NGC3718, since the observa¬ 
tions appear to be compatible with key as pects of the current the- 
oretical view around SMBH recoils (e.g.lKomossa 2012; Sijacki 
|etal.| 201 l| ). 

24 This velocity shift represents the projected kick velocity, so it may 
only serve as a lower limit of the true kick velocity and this should be 
accounted for, in all subsequent calculations. 
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7. Discussion 

Summarizing the results, the observations suggest that 
NGC3718 appears to have, a spatial offset which, potentially co¬ 
incides, with a radio emission offset, a relative redshift between 
its broad and narrow Ha lines and the presence of a core in its 
surface brightness profile, all compatible with the presence of an 
SMBH recoil. Moreover, the strong indications that NGC3718 is 
a merger remnant including, spiral-like photometry mixed with 
elliptical-like kinematics and properties that belong to both core 
and coreless ellipticals, fulfill the prerequisite for an SMBH re¬ 
coil, in the sense that merging is the current theoretical mecha¬ 
nism for the formation of an SMBH binary that can lead to an 
SMBH recoil. Some of these indications, however, could also be 
viewed to apparently contradict each other. Thus, a very impor¬ 
tant question is the following: “Is there a theoretical framework 
which could incorporate all the, potentially contradictory, char¬ 
acteristics, into one physically consistent scheme?”. 


7.1. Formation of cores and galaxy evolution 

Attempting to give a decent answer to this question, requires us 
to briefly review what we know, up to now, about the two rele¬ 
vant sub-classes of the elliptical galaxies family. In recent years, 
it has become widely accepted that, within the A CDM cosmol- 
ogy, structures grow hierarchically in the Cosmos (e.g. |White| 
& Rees |1978| . Within this picture, classical bulges and ellipti- 
cal galaxies are formed in major mergers, resulting in the com¬ 
mon properties observed in these objects. Moreover, elliptical 
galaxies define a trend of increasing dissipation with decreasing 
mass ( |Kormendy|2009| ). This means that later-type, fainter ellip¬ 
ticals, being “1st stage” merger products of gas-rich spiral galax- 
ies, have mo re gas to dissipate through “wet” merging (Hopkins 
et al. |2009a| , whereas earlier type giant ellipticals, being prod¬ 
ucts of the subsequent merging of gas-poor fainter ellipticals, 
were formed in “dry” mergers ( Hopkins et al.|2009b| ). 

An immediate effect of “wet” mergers is that older stars 
(formed in the progenitors before merging), relax violently at 
larger radii (giving rise to a Sersic law profile), while new stars 
are formed mainly in a central intense starburst due to dissipa¬ 
tion. As the cold gas experiences tidal torques due to the violence 
of a major merger, is channeled towards the central region, trig¬ 
gering SMBH growth, SMBH feedback and an intense nuclear 
starburst ( [Hopkins et al.|2008| . The latter forms the “extra light” 
component in these objects, which is observationally recogniz¬ 
able as a break towards a steeper inner surface brightness profile, 
when compared to an inward extrapolation of the outer profile’s 
Sersic fit (e.g. |Hopkins et al.||2009a|b[ [Kormendy & Ho[|2013| ). 
During this process the SMBH accretes up to the point that the 
energy feedback output from the AGN is large enough (~ 1 % of 
a, near the Eddington limit, AGN’s energy output is coupled to 
the infalling gas) to blow away the residual gas, stopping the BH 
growth and quenching the star formation in the central region 
(e.g. |Ostriker & Ciotti|2005| ). SMBH recoils are still expected to 
occur, since merging is a prerequisite of the process of ellipti¬ 
cal galaxies formation. The presence of gas, however, helps the 
SMBH recoil to return to the center more quickly (e.g. Blecha| 
|et al.|20lH|Sijacki et al.|20lT] ), stopping any core scouring and 
replacing the “missing light” with newly formed stars in the cen¬ 
tral starburst ( [Kormendy & Ho|2013| ), which in turn, consumes 
the larger fraction of the available gas (Hopkins et al.| |2Q08). 
During this process, the AGN operates in quasar mode. 

However, “dry” mergers, as a product of merging of fainter 
gas-poor “extra light” ellipticals, lack the intense central star- 


burst. The old violently relaxed stellar populations of the pro¬ 
genitors will again violently relax at larger radii in the newly 
formed system, whereas the “extra light” components will end 
up near the center and be preserved, in the form of a m ore 
compact stellar distribution in the central region ( [Hopkins et al. 
2009bj ). Although tidal forces will channel the (little) remaining 


gas towards the center, the newly formed giant elliptical is mas¬ 
sive enough to hold onto enough amounts of X-ray emitting gas, 
which ensures that any cold gas reaching the nuclear region will 
be heated up, preventing any star formation (Nipoti & Binney 


2007). The AGN in this case operates in the, so called, mainte¬ 


nance mode, which further helps the gas to remain hot, making 
the nuclear star formation even more difficult and , eventually, 
keeping the merger “dry” ( [Kormendy & Ho||2013| . Due to the 
absence of gas, SMBH recoils, are able to “stay on duty” for 
a longer time, acting mainly on a spacial range of the order of 
the preserved “extra light” fossil in the newly formed remnant, 
scouring the core and creating what appears as a “missing light” 
or, in other words, a break towards a shallower inner surface 
brightness profile, when compared to an inward extrapolation of 
the outer profile’s Sersic fit. 

7.2. Formation of polar ring galaxies 

As we extensively discuss in section HU the observations sup¬ 
port NGC3718’s classification as a polar ring galaxy, suggested 
by |Sparke et al.| (|2009). But how are polar ring galaxies formed? 

There are two candidate formation mechanisms: The merg¬ 
ing scenario, proposed by Bekki ( 1997| 1998| ), suggesting for¬ 
mation through a head-on collision between two orthogonally 
placed spiral galaxies, the intruder and the victim (left panel 
of their figure 1). After the merging the intruder becomes the 
host galaxy (often an SO), whereas the victim becomes the po¬ 
lar ring. The second is the accretion scenario. Supported by a 
number of authors (e.g. Reshetnikov & Sotnikova|[1997| ), it ex¬ 
plains polar ring galaxy formation through gas accretion from a 
donor galaxy, as a result of tidal gravitational interactions due 
to a close encounter. Subsequent merging of the two galaxies is 
not necessary for the formation of a ring, but it is not physically 
prohibited. 

[Bournaud & Combes] ( |2003[ ) simulated both of these ideas, 
using the same numerical model. They find that both mecha¬ 
nisms successfully reproduce many of the observed properties of 
polar ring galaxies such as, a) early-type host galaxies, b) stable, 
nearly polar, ring structures (some of their runs are evolved for 
8-10 Gyr with a typical formation time for a polar ring galaxy 
of ~ 3 Gyr) and c) warps/spiral arm like structures when these 
objects are viewed from specific lines of sight (edge-on/face-on 
for the ring/host galaxy respectively). All of these properties are 
observed to be pr esent in NGC3718, by a number of indepen- 
dent studies (e.g. Pott et al.l|2004[|Krips et al.|2005||Sparke et al. 
[2009] ). 

A few properties, however, differentiate the two proposed 
scenarios. Firstly, in the merging scenario the host galaxy has 
to be substantially free of gas (which is the case of NGC3718 
according to [Sparke et al. ( |2009| ), whereas in the case of the 
accretion scenario the host galaxy can also be a gas-rich ob¬ 
ject. Secondly, the merging scenario is expected to form a faint, 
spherical, diffuse stellar halo, consisting of the stars of the vic¬ 
tim galaxy and surrounding both the host galaxy and the polar 
ring. Such a faint diffuse stellar component can be seen in fig¬ 
ure [T] of NGC3718, nearly engulfing the edge-on ring and the 
host galaxy. However, in the accretion scenario, instead of such 
a component, a donor galaxy has to be identified in a relatively 
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close proximity, which in this case could be the close neighbor¬ 
ing galaxy NGC3729 or a similar galaxy that has already been 
cannibalized by NGC3718. Lastly, the merging scenario appears 
to be more consistent, in (inefficiently) transferring gas towards 
the central region (~ 10 - 25% in most of their runs), than the 
accretion scenario. [Boumaud & Combesj ( |2003| ) quote that, the 
gas infall is not as large and as general in the accretion scenario 
as it is in the merging scenario. They conclude that the accretion 
scenario is the most likely formation mechanism, at least, for the 
majority of the cases without, however, ruling out the merging 
one based on some physical explanation, but rather through sta¬ 
tistical argumentation, leaving space for it to actually occur in 
Nature. 


To the extent that a subsequent merging, after a possible ac¬ 
cretion of gas into polar orbit, is not a common process, our 
analysis tends to favor the merging scenario for the formation 
of NGC3718. The reason for this is mainly the observational ev¬ 
idence we present throughout this work for the presence of an 
SMBH recoil, which requires a merger to occur. Moreover, the 
low star formation activity (also seen by |Sparke et al. ( 2009 }) and 
the detection of rap idly rotating molecu lar gas within the central 
~ 700 pc seen by [Krips et aL| ( ]2005| (through their position- 
velocity diagrams), also favor the merging scenario, since the 
gas inflow predicted by this appears to be consistent with these 
observations. Finally, if the dust component traces the kinemat¬ 
ics of the gas (a reasonable assumption since, gas and dust are 
mixed in the apparent polar ring), then the contribution from hot 
dust to the K s band emission that we detect (see section |4.1| ), 
serves as an additional indirect tracer for the aforementioned gas 
inflow. 


The accretion scenario can also, in principal, explain these 
phenomena. Provided that a subsequent merging is a common 
outcome of this scenario, then this could be viewed as a merg¬ 
ing event of different orbital configurations. So, when we refer 
to a “merging event”, we do mean either of the aforementioned 
polar ring galaxy formation scenarios, as long as the merging of 
SMBHs is “guarantied”, in order to explain the indicated SMBH 
recoil and the observed core. 


that can flow towards the center, triggering the necessary star 
formation that could “fill the gap” created by SMBH scouring. 

Additionally, the behavior of the A4 parameter (figure |T5] ), 
indicative of an inside-out transformation from a disky to a boxy 
object, fits naturally to this framework. Boxyness is present in 
objects that lack the intense central starburst, either because 
of limited amounts of nuclear cold gas and/or star formation 
quenching due to AGN feedback ( [Kormendy & Ho [2013 ) . Such 
an environment seems ideal for an SMBH recoil to excavate a 
core. This process appears to be connected with the isophote 
contours shape behavior in the case of NGC3718, since the tran¬ 
sition from boxy to disky, occurs roughly in the spatial scales of 
the observed core. 


In this framework, the presence of energetic X-ray emission 
(14 - 195 Kev) reported by |Tueller et al.[ ( [20T0] ), also fills a gap. 
Provided that it is associated with AGN feedback, it makes the 
star formation quenching mechanism applicable in this case. The 
very good agreement between the velocity dispersion measured 



Any apparent contradictions originate from the fact that, all 
these processes do not take place in a “red and dead” giant ellip¬ 
tical (for which they were originally proposed) but, circumstan¬ 
tially, due to the limited gas inflow attributed to the nature of the 
dynamics of this (rare) merger case, in a gas-rich system. 

The last of NGC3718’s properties that we have to incorpo¬ 
rate into this scheme, is its spectral classification. NGC3718 is 
classified as a LINER (low-ionization nuclear emission line re¬ 
gion) galaxy ( Ho et al. 1997b) ) and it is also part of the NUGA 
sources, a survey aimed at the stu dy of nearby low-luminosity 
active galactic nuclei (LLAGNs) (Garcia-Burillo et al. |2003| ). 
This is also consistent with the inefficient transport of gas to¬ 
wards the center, in the sense that it affects the accretion rate 
by reducing the available gas reservoir. The latter is implied by 
the fact that, NGC3718 appears to be a sub -Eddington system 
with ~ 10 -4 - 10 -5 (Krips et al. 2007), indicating that the 
accretion rate onto the central engine is fir from being efficient. 


7.3. Putting the pieces together 


7.4. On the SMBH recoil 


The inefficient transport of gas towards the nucleus, however, 
might be the key in providing a robust physical explanation for 
the core that we see in the surface brightness profile of figure [T4| 
A working scheme could be the following: 

If NGC3718 is viewed as a gas-richp] merger remnant, one 
would expect to see signs of the intense nuclear starburst in 
the form of an “extra l ight” component, as describ ed by, e.g. 
Kormendy & Ho| ( [20T3] ); |Hopkins et al.| (|2008 12009a), instead of 
the observed core. The suggested absence of sufficient amounts 
of nuclear gas, however, lead naturally to the limited star forma¬ 
tion observed. This, combined with the fact that in the absence of 
sufficient amounts of gas, the time scale for which an SMBH re¬ 
coil “stays on duty” is prolonged ( [Sijacki et al.|2011| ), provides 
a physically consistent framework for the formation of the ob¬ 
served core. Put more simply, the gas exists in the galaxy but it 
is “locked up^_on (the stable) polar orbit ( Bournaud & Combes 
2003) |Bekki||1997 |1998| , a fact that reduces the available gas 


25 Sparke et al. ([2009) measures M m ~ 8 x 10 9 M 0 , a value ~ 2 
times larger than the Milky Way, while Krips et al. (2005) measures 
~ 2.4 x 10 8 M o of molecular gas through CO observations, values that 
can characterize NGC3718 as a gas-rich galaxy. 


The last question within the scope of this study, is, perhaps, the 
most prevalent one, albeit a very difficult one to answer: “Is this 
indeed a case of a true SMBH recoil?”. An attempt to quali¬ 
tatively answer this question requires us to further discuss the 
compatibility of an SMBH recoil with other observational find¬ 
ings for NGC3718. 

The near-the-center position of the offset red blob, allows 
space for two different scenarios: It could either represent the 
beginning of an SMBH displacement, or, the “capture” of an 
ongoing SMBH recoil in one of its pericentric passages. 

At first, if the age of ~ 2 - 3 Gyr of NGC3718’s gas disc, 
suggested by [Sparke et al. ( 2009| ), traces roughly the age of the 
merger, then it is in excellent agreement with the results of the 
simulations of [Boumaud & Combes ( 2003[ ), where a stable polar 
ring and a warp is formed within ~ 2 - 3 Gyr, from mergers with 
mass ratios of ~ (1-4) : l. |Blecha et al.| ( |2Ql7] ) simulated recoils 
in galaxy mergers and they find that, for mergers with mass ratios 
~ (1 - 2) : 1 and gas content of 10 - 30%, the coalescence of 
the SMBHs often occurs at ~ 2 Gyr. For such systems and for 
kick velocities of v^ ic k < 0.1v esc , the trajectories of the kicked 
SMBHs are often confined within ~ 1 kpc with setback times of 
the order of ~ 1 Gyr. So, time wise, it seems plausible that either 
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a new or an ongoing SMBH recoil exists in NGC3718 and that 
it would be observable. 


Moreover, assuming that our current view for the formation 
of cores is attributed mainly to SMBH binaries and the subse¬ 
quent recoils (e.g. Kormendy & Hoj20T3 ), then the rather large 
Mdef ~ 19 Mbh ( see section |5.2|), needs to be addressed. Such 


a large M de f could be partly attributed to accumulated errors 
from i.e. uncertainties in the fit, uncertainties in the ^ relation 
and/or due to the large intrinsic scatter of the scaling relations. 
Such large mass deficits, however, could also imply that pre¬ 
existing co res in the progenitor galaxie s have been inherited to 
NGC3718. |Gualandris & Merritt ( |2008| ) simulated kicks and the 
mass deficits that these can pr oduce. The y find that a combina¬ 
tion of a pre coalesced binary (Merrit tp006| ) followed by a sub¬ 
sequent recoil, can excavate as much as ~ 5 M BH per merger. The 
fact that the light profiles are (more or l ess) p reserved in dissipa¬ 
tionless mergers (e.g. [Hopkins et al.|2009b| ), is what leads sev¬ 
eral authors to the conclusion that, large mass deficits could be 
indicative of different core excavation events following sequen¬ 
tial mergers. A large M de f therefore, seems possible to co-exist, 
with either a new or an ongoing SMBH recoil. 

NGC3718’s radio emission, however, may help us draw a 
line between a new and an ongoing SMBH recoil. To the extent 
that NGC3718 is indeed an X-shape source, it could be the host 
of a recently recoiled SMBH that has undergone a reorientation 
of its jet (Liu et al. 2012 ). In this case, the inheritance of a core 
of the observed size (~ 200 pc across) is considered necessary, 
since a just recently recoiled SMBH would not have the time to 
excavate such a large (both spatially and in terms of M de f ) core 
(Merri tt|2006| ). If, on the contrary, the secondary extensions on 
the radio map are not attributed to a second jet-fossil, then, an 
ongoing SMBH recoil which accretes during one of its pericen¬ 
tric passages ( [Blecha etaL 201 1| ), could produce a single bipolar 
jet structure. In this case, the inheritance of a core of the ob¬ 
served size is not necessary (but, partly, not prohibited either). 
So, in principal, a core could be excavated in “one go”, provided 
that the M de f is considerably overestimated due to the aforemen¬ 
tioned uncertainties. 


To conclude this work, we have presented much indepen¬ 
dent evidence, which, along with our measurements on the 
SUBARU and e-Merlin data, suggest that it is reasonable to treat 
NGC3718 as a very good candidate host galaxy of an SMBH re¬ 
coil. Therefore, an interesting question arises: “Is the presence 
of SMBH recoils in mergers, a standard process in galaxy evo¬ 
lution?”. If, as the observations suggest, NGC3718 is indeed 
a polar ring galaxy, and polar ring galaxies make up ~ 5% 
of the lenticular and early type galaxies in the local universe 
( [Whitmore et al.j jl 990ft , then, the detection of an SMBH recoil 
in such a rare object may, by itself, suggest an affirmative an¬ 
swer to the above question. In any case, this is as far as our data 
allow us to go in interpretation terms. Our limited knowledge 
on the exact behavior of recoiled SMBHs, our incomplete un¬ 
derstanding of the precise mechanism(s) that are responsible for 
the formation of cores and the limited resolution of the current 
observations, allow space only for general qualitative, and there¬ 
fore highly speculative, interpretations. Certainly though, bet¬ 
ter quality future observations combined with higher resolution 
simulations, are necessary, in order to answer the above ques¬ 
tions in a more robust, quantitative way. In all, NGC3718 turns 
out to be an object full of wonderful surprises with, possibly, a 
lot more that are yet to be revealed. 


8. Summary 

We have extensively studied NGC3718, using NIR SUBARU 
and 6cm e-Merlin radio data as well as previously published re¬ 
sults. Our findings are summarized as follows: 

1. Our NIR color maps do not show any signs of a large scale 
constant color gradient, indicating that our view towards the 
nucleus of NGC3718 is unaffected by extinction. 

2. An offset red blob is detected in our NIR color maps, being 
displaced by ~ 4.25 pc from the center of the underlying 
stellar bulge. The radio emission appears to originate from a 
position closer to the offset red blob than to the center of the 
stellar bulge. 

3. A light decomposition reveals a contribution (up to ~ 50%) 
from a hot (up to ~ 1000K) dust component to the incoming 
light from the central ~ 0.5 arcsec, indicative of the presence 
of gas in that region. 

4. An extended elongated structure ~ 1 arcsec across, is de¬ 
tected in the e-Merlin radio map, probably indicative of a 
small scale bipolar jet. A second smaller one ~ 0.6 arc¬ 
sec across, lies almost perpendicular to the first one, sug¬ 
gesting a possible X-shape radio source. However, these ex¬ 
tended structures may also contain contribution from SNR 
(Supernova Remnant) related non-thermal emission. 

5. A disagreement between the M B h - L and M B h - cr scal¬ 
ing relations when NGC3718’s measured values are used and 
the shape of the surface brightness profile, reveal an amount 
of missing light in the form of a scoured core of ~ 200 
pc across. This light “deficit” translated to mass, reflects an 
M de f ~ 19 Mbh- 

6 . NGC3718 shows mixed characteristics in two ways: On 
one hand, it has spiral-like photometry (seen in the sur¬ 
face brightness profile as a dual disk component) combined 
with elliptical-like kinematics (seen in the surface bright¬ 
ness profile as a bulge component with Sersic index n > 4). 
Furthermore, NGC3718 shows characteristics that belong to 
both, fainter, more rotationally supported, disky, coreless el¬ 
lipticals, as well as to giant, less rotationally supported, boxy, 
core ellipticals. These, strongly imply that NGC3718 is the 
result of a merger. 

7. Finally, NGC3718 appears to be the host of an SMBH re¬ 
coil. Evidence for that include an offset NIR red blob, with 
a, potentially coincident, offset radio emission, the presence 
of a core in the surface brightness profile, a relative redshift 
between the broad and narrow Ha emission lines and the in¬ 
dication for the presence of a Hyper Compact Stellar System 
surrounding the offset NIR red blob. 
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Appendix A: On the influence of the dust lane 

Appendix [A| provides additional evidence to the ones discussed 
in section ^regarding the potential influence of the dust lane on 
our view towards NGC3178’s nucleus. 

For this reason we present figure |A.1| where the 2MASS 
K s band image of NGC3718 is displayed. What can be noted 
in figure [ATT] is the fact that the 2MASS Ks band contours (in 
red) appear to be, almost entirely, unaffected by the presence of 
the dust lane. They appear to be significantly rounder and bet¬ 
ter defined, especially when compared to the SDSS r band con¬ 
tours (in green), which are severely deformed by the presence 
of foreground dust in the region. Additionally, our SUBARU K s 
band contours (in yellow) appear to be more of a smaller scale 
inward extrapolation of the 2MASS Ks band contours (as ex¬ 
pected), rather than of the visible light, heavily, distorted con¬ 
tours. Moreover, the innermost - 4 arcsec of our SUBARU Ks 
band, are located ~1.5 arcsec away from the region where the 
contours deformation in the optical becomes catastrophic, which 
naturally, coincides approximately with the projected beginning 
of the dust lane (white dashed box). 

We consider as safe, therefore, to treat the central region of 
NGC3718 in the NIR, as being unaffected by the presence of the 
dust lane, when we attempt to interpret our NIR color maps. 


Appendix B: On the z-J bands alignment 

Appendix [B] focuses on the validity of the coordinates calibra 
tion method used on the SUBARU data, described in section[3U 


For this reason we present figure B.l where the SDSS z 
band image of NGC3718 is displayed. In this image, the black 
dashed ellipse has a major axis of - 4 arcsec, while the blue 
dashed line indicates that this is, approximately, the radial dis¬ 
tance where the SDSS z band contours begin to deviate from 
being symmetric with respect to the peak flux positions. What 
can be seen, is that the SDSS z band contours are well defined, 
relatively round and, more importantly, symmetric within this 
region. Especially when compared to the SDSS r band contours 
(green contours in figure [AT] ), they seem to be almost entirely 
unaffected by the presence of dust within the region of interest 
and we consider them to be suitable for use, for an accurate es¬ 
timation of the photocenter of NGC3718 in this band. Moreover 
our J band contours appear to be similar to the SDSS z band 
contours, i.e. in terms of similar roundness. The z band elliptic- 
ities range between -0.13- 0.15, while in J band they range 
between - 0.09 - 0.12, in the radial interval - 0.5 - 1.5 arcsec, 
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Figure A.l. 2MASS K s band image of NGC3718. Over-plotted 
can be seen: 1) 2MASS K s band contours (red), 2) SDSS r band 
(~ 0.6 /d m) contours (green), 3) SUBARU K s band contours 
(yellow) and 4) distance (blue dashed line length ~ 3.5 arcsec) 
of the optical and NIR peak flux positions to the approximate 
beginning of the dust lane (outlined with the white dashed box). 


indicating that, at least to a large extent, both bands appear to be 
similarly unaffected by the presence of foreground dust. 

We consider, therefore, that the aforementioned symmetry of 
the SDSS z band contours in the central ~ 4 arcsec of NGC3718, 
as well as their similarity to the J band contours, sa fely allow us 
to apply the centering method described in section |3.1| The de¬ 
rived photocenters are used to align the SDSS z band with our 
SUBARU J,H and K s bands, in order to astrometrically cali¬ 
brate them. 


H-Ks 



R.A. (11:32:34.sss) 


Figure C.l. Same as figure [l3| but with the radio map position 
dictated by the old phase reference source position. Thisj_s the 
equivalent position to the one reported by Krips et al. (2007]), 
namely R.A.: 11:32:34.8534 ± 0.0005, Dec.: +53:04:04.523 ± 


0.004. 



Figure B.l. SDSS z band image of NGC3718. Over-plotted can 
be seen: 1) The SDSS z band contours (yellow) and 2) the 
SUBARU J band contours (red). The black dashed ellipse and 
the blue dashed line show the symmetry of the SDSS z band 
contours within the central ~ 4 arcsec of NGC3718. 


Appendix C: On the e-Merlin radio map position 
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